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1 

METHOD AND APPARATUS FOR DETERMINING MUSICAL NOTES 

FROM SOUNDS 

Field of the Invention 

5 

The present invention relates to determining musical notes from sounds, such 
as humming or singing. In particular it relates to converting such sounds into 
notes and recognising them for the purpose of music retrieval. It also relates 
to the component means and processes. 

10 

Background Art 

Multimedia content is an increasingly popular resource, supported by a 
surging market for personal digital music devices, an increase of bandwidth to 

15 the home and the emergence of 3G wireless devices. There is an increasing 
need for an effective searching mechanism for multimedia content. Though 
many systems exist for content-based retrieval of Images, few mechanisms 
are available to retrieve the audio portion of multimedia content. One 
possibility for such mechanisms is retrieval by humming, whereby a user 

20 * searches by humming melodies of a desired musical piece into a system. 
This Incorporates a melody transcription technique. 

Figure 1 shows a flowchart for a known system of humming recognition. The 
melody transcription technique consists of a silence discriminator 101, pitch 
25 detector 102 and note extractor 103. It is assumed that each note will be 
separated by a reasonable amount of silence. This reduces the problem of 
segmentation to a silence detection problem. 

In U.S. Pat. No. 6,188,010 a FFT (Fast Fourier Transfonn) algorithm is used 
30 to analyse sound by obtaining the frequency spectrum information from 
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waveform data. The frequency of the voice is obtained and finally a music 
note that has the nearest pitch is selected. 



In U.S. Pat. No, 5,874,686 an autoconrelation-based method is used to detect 
5 the pitch of each note. In order to improve the perfonmance and robustness 
of the pitch-tracking algorithm, a cubic-spline wavelet transfomi or other 
suitable wavelet transform is used. 



In U.S. Pat No. 6,121.530 the onset time of the voiced sound is divided off as 
10 an onset time of each note, a time difference with an onset time of the next 
note is determined as the span of the note and the maximum value among 
the fundamental frequencies of each note contained during its span is defined 
as the highest pitch values. 

IS Automatic melody transcription is the extraction of an acceptable musical 
description from humming. Typical humming signal consists of a sequence of 
audible waveforms interspersed with silence. However, there is difficulty in 
defining the boundary of each note in an acoustic wave and there Is also 
considerable controversy over exactly what pitch is. Sound recognition 

20 involves using approximations. Where boundaries between notes are clear 
and pitch is constant, the prior art can produce reasonable results. However, 
that is not necessarily so where each audible waveform may contain several 
notes and pitch is not necessarily maintained, as happens with real people 
humming. A hummer's inability to maintain a pitch often results in pitch 

25 changes within a single note, which may be subsequently misinterpreted as a 
note change. On the other hand, if a hummer does not pause adequately 
when humming a string of the same notes, the transcription system might 
interpret It as one note. The task becomes increasingly difficult in the 
presence of expressive variations and the physical limitation of the human 

30 vocal system. 
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Object and Summary of the Invention 

It is therefore an aim of the present invention to provide an improved system 
5 for recognising hummed tunes or the like and to provide component 
processes and apparatus that can be used in such a venture. 

According to a first aspect of the invention, there is provided a method fjor use 
In transcribing a musical sound signal to musical notes, comprising the steps 
10 of: 

producing note markers, indicative of the beginnings and ends of notes 
in said sound signal; and 

detecting the pitch values of notes marked by said note markers. 

15 Preferably this method further comprises detecting portions of said sound 
signal that can be deemed to be silences. 

This method may also further comprise the step of extracting notes from said 
pitch values to create note descriptors. 

20 

According to a second aspect of the invention, there is provided a method for 

detecting portions of a musical sound signal that can be deemed to be 

silences, comprising the steps of: 

dividing said sound signal into at least one group of blocks; 
25 deriving short-time energy values of said blocks in a group; 

deriving a threshold value based on said short-time energy values; and 
using said threshold value to classify blocks of said group as silent or 

otherwise. 
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According to a third aspect of the invention, there is provided a method of 
producing note markers, indicative of the beginnings and endings of notes in 
a musical sound signal, comprising the steps of: 

extracting an envelope of said sound signal; 
5 differentiating said envelope to compute a gradient function; and 

extracting note markers from said gradient function, indicative of the 
beginnings and ends of notes in said sound signal. 



The process of envelope extraction may comprise the steps of: 
10 performing full-wave rectification on said sound signal; and 

low-pass filtering the output of the full-wave rectification. 



The process of differentiation may comprise the steps of: 
determining the gradient of said envelope; and 
15 low-pass filtering said gradient 



The process of note markers extraction may comprises the steps of: 
removing small gradients from said gradient function; 
extracting turning points of the attack and decay of remaining 
20 gradients; 

removing unwanted attacks and decays; and 

registering remaining attacks and decays as said note mariners. 

According to a fourth aspect of the invention, there is provided a method for 
25 detecting the pitch values of notes in a musical sound signal, comprising the 
steps of: 

isolating notes in the sound signal; 

dividing said notes into one or more groups of blocks; 

deriving pitch values of said blocks; and 
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deriving the pitch values of said notes by means of clustering on said 
pitch values of said blocks. 

This process of Isolating notes may use note mariners to do so. 

5 

One or more of the above aspects may be combined. 

According to a fifth aspect of the invention, there is provided a method of 
identifying pieces of music, comprising the steps of: 
10 receiving a musical sound signal imitative of a piece of music; 

transcribing said musical sound signal to a series of musical notes and 
timings using the method of the first aspect above; 
comparing said series of musical notes and timings with series of notes 
and timings of pieces of music in a database; and 
15 identifying the piece of music deemed most similar by this comparison. 

Following this, the identified piece of music may then be retrieved. 

The invention is not limited to human use. It may be useful in conducting 
20 experiments with animals. Moreover, it is not limited to humming, but could 
be used with whistling, singing or other noise production. 

The Invention also provides apparatus operable according to the above 
methods and apparatus corresponding to the above methods. 

25 

This method and apparatus extract symbolic high-level musical structure 
resembling that of a music score. Humming or the like is converted with this 
invention into a sequence of notes that represent the melody that the user 
(usually human, but potentially animal) is trying to express. These retrieved 
30 notes each contain information such as a pitch, the start time and duration 
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and the series contains the relative order of each note. A possible application 
of the invention is a music retrieval system whereby humming fomis the query 
to some search engine. Music retrieval via query-by-humming can be applied 
to different applications such as PC, cellular phone, portable jukebox, music 
S kiosk and car jukebox. 

Brief Description of Drawing 

10 The present Invention is now further described by way of non-limitative 
example, with reference to the accompanying drawings, in which:- 

Figure 1 is a flowchart of a prior art melody transcription technique; 

15 Figure 2 is a schematic block diagram of an embodiment of the present 
invention; 

Figure 3 is a flowchart of a melody transcription technique used in the 
embodiment of Figure 2; 

20 

Figure 4 is a flowchart of operation of a silence discriminator used in the 
embodiment of Figure 2; 

Figure 5A is a flowchart of gradient-based segmentation used in the 
25 embodiment of Figure 2; 

Figure 5B is an illustration of a typical humming wavefomi; 

Figure 5C is an illustration of the output of the envelope detector, with the 
30 waveform of Figure 5B as the input; 
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Figure 5D is an illustration of the output of the differentiator, with the 
wavefonn of Figure 5G as the input; 

Figure 5E is an illustration of the note markers produced by the note markers 
5 extractor, with the waveform of Figure 5D as the input. 

Figure 6 is a flowchart of operation of an envelope detector used in the 
embodiment of Figure 2; 

10 Figure 7 is a flowchart of operation of a differentiator used in the embodiment 
of Figure 2; 

Figure 8 is a schematic illustration of the criteria for selection of legitimate 
attack and decay; 

15 

Figure 9 is a flowchart of operation of a note marker extractor used In the 
embodiment of Figure 2; 

Figure 10 is a flowchart of a gradient threshold function used in the 
20 embodiment of Figure 2; 

Figure 11 Is a flowchart of operation of an edge detector used in the 
embodiment of Figure 2; 

25 Figure 12 is a flowchart of operation of a pitch detector used in the 
embodiment of Figure 2; and 



30 



Figure 13 is a flowchart of operation of a dominant pitch detector used in the 
embodiment of Figure 2. 
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Specific Description 

A robust melody transcription system is proposed to serve as an ensemble of 
solutions to solve the problem of transcribing humming signal to note 
5 descriptors. A melody technique is used to produce note descriptors. This 
Information is used by a feature extractor to obtain features to be used in a 
search engine. 

Figure 2 is a schematic block diagram of an embodiment of the present 
10 invention. A digitised humming input signal S200 from a PC, cell phone, 
portable Jukebox, music kiosk or the like, is input into a melody transcription 
device 2. There it is input In parallel Into a pitch detector 202, a silence 
discriminator 204 and a gradient based segmentation unit 206, where it first 
goes into an envelope detector 208. The envelope detector 208 produces an 
15 envelope signal S210 from the humming signal, which is input into a 
differentiating circuit 212. Another input into this is a silence mariner signal 
S214 from the silence discriminator 204. The output from the differentiating 
circuit 212 is a gradient function signal S216, which is input Into a note marker 
extractor 218, which also receives the silence mari<er signal S214 from the 
20 silence discriminator 204. The note mariner extractor 218 outputs a note 
mariner signal S220, which, together with the silence marker signal S214 and 
humming input signal S200, is input Into the pitch detector 202. The gradient 
based segmentation unit 206 is made up of the envelope detector 208, the 
differentiating circuit 212 and the note marker extractor 218. 

25 

Using the three inputs, the pitch detector 202 produces a pitch value signal 
S222, from which a note extractor circuit 224 produces a note descriptor 
signal S226. This then is output from the melody transcription device 2. In 
this example, a feature extraction circuit 228 produces a feature signal S230, 
30 from the note descriptor signal S226. An MPEG-7 descriptor generator 232 
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uses this to produce a feature descriptor signal S234. which Is fed to a search 
engine 236. Searching using a music database 238 gives a search result 
S240. 

5 The silence discriminator 204 illustrated in Figure 2 is employed to Isolate the 
audible portion of the input humming signal S200 from the silence. The pitch 
detector 202 is used to compute the pitch of the humming input S200. The 
structure of the audible waveform is complex but the present invention uses 
detection of an attack and decay pair indicating the existence of a note. Thus 

10 the envelope detector 208 is employed to remove the complex structure of 
the audible wavefomn. The differentiator 212 computes the gradient of the 
envelope S210. Another difficulty is the ambiguous nature of the attack and 
decay pair that symbolises the existence of a note. Unlike musical 
instruments, people cannot transit to the next note with a boundary that Is well 

15 defined. The problem Is compounded by the fact that the volume may 
change due to expression or failure by the hummer to maintain the volume. 
The volume change might create a false attack and decay within the duration 
of a particular note. The note marker extractor 218 is therefore used to 
remove all the ^Ise attacks and decays. The legitimate attack and decay 

20 pairs left are used as note makers that mark the start and end of a note. With 
the knowledge of the location of each note, the pitch detector 202 computes 
the pitch of each note. Finally, the note extractor 224 is employed to map the 
pitch values and note markers to produce note descriptors. A note descriptor 
contains infomnation such as pitch, start time and Interval of a particular note. 

25 

In this preferred embodiment, the melody transcription system comprises two 
distinct steps: segmentation and pitch detection. The segmentation step 
searches the digital signal S200 to find the start and duration of all notes that 
the hummer tries to express. The silence discriminator 204 isolates the 
30 voiced portions. This infonnation has been used in the prior art to segment 
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the digital signal. This Is only feasible if a hummer inserts a certain amount of 
silence between each note. Most inexperienced hummers have difficulties 
Inserting silence between notes. In this invention, a gradient-based 
segmentation method is employed to search for notes within the voiced 
5 portions, thus not relying so much on silence discrimination. 

The humming signal is similar to an amplitude modulated (AM) signal where 
the volume is modulated by the pitch frequency. The pitch signal Is not useful 
in this case, which Is removed to extract the envelope. The envelope shows 

10 some interesting properties of a typical humming signal. The envelope 
increases sharply from silence to a stable level. The stable level is 
maintained for a while before it drops back sharply to silence again. Thus the 
existence of an attack, followed by a steady level and a decay of a note, is 
evidence of the existence of a note. The gradient-based segmentation is 

15 derived from these unique properties to extract the note markers. 

These note markers are used in this invention to enhance the performance of 
the pitch detector 202. The approach is to exploit the fact that the pitch within 
each pair of start and end note makers is supposed to be constant. The 

20 signal of each note is divided Into blocks of equal length. The signal in each 
block is assumed to be stationary and the pitch (frequency) is detected by 
autocon'elatlon. In an ideal case, these values are identical. However, the 
autocorrelation pitch detector 202 is sensitive to harmonics that cause errors 
in the detection of pitch. Furthemiore, hummers frequently fail to maintain the 

25 pitch within the duration of a particular note. A k-mean clustering algorithm is 
selected in this Invention to find the prominent pitch value. 

Music retrieval by humming is perceived as an excellent complement to tactile 
interfaces on handheld devices, such as mobile phones and portable 
30 jukeboxes. This Invention can also be employed in a ring-tone retrieval 
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system whereby a user can download the desired ring-tone by humming to a 
mobile device. 



Thus, in this embodiment, a user hums a tune into a microphone attached to 
5 a PC, cell phone, portable jukebox, music l<iosl< or the like, where the input 
sound is converted into a digital signal and transmitted as part of a query. 
The query is sent to a search engine. Melody transcription and feature 
extraction modules in the search engine extract relevant features. At the 
same time, the search engine requests MPEG-7 compliant music metadata 

10 from music metadata servers on its list. The search proceeds to match the 
music metadata with the features extracted from the humming query. The 
result is sent back to the user, with an indication of the degree of match (in 
the form of a score) and the location of the song(s). The user can then 
activate a link provided by the search engine to download or stream the song 

15 from the relevant music collection server - possibly for a price. The MPEG-7 
descriptor generator is optional and depends on the application scenario. 

Such a mechanism entails a robust melody transcription subsystem, which 
extracts symbolic high-level musical structure resembling that on a music 
20 score. Thus the humming must be converted into a sequence of notes that 
represent the melody that the user tries to express. The notes contain 
information such as the pitch, the start time and the duration of the respective 
notes. Thus it requires two distinct steps: the segmentation of the acoustic 
wave and detection of the pitch of each segment. 

25 

In the prior art shown in Figure 1, the melody transcription technique consists 
of a silence discriminator, pitch detector and note exb-actor. Figure 3 is a 
similar flowchart, showing the components of the present invention. Once 
again, there is a silence discriminator step 301 and a pitch detector step 304, 
30 which leads to a note extractor step 305. However, in this invention, an 
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additional step is introduced into tlie conventional technique In the fonm of an 
•advanced mode' option step 302, following on from silence discriminator step 
301. The selection of the advanced mode activates the gradient-based 
segmentation step 303. This step Is made up of the processes conducted In 
5 the gradient based segmentation unit 206 of Figure 2. Thus the process 303 
searches for note markers within each voiced wavefomi. Note markers found 
are processed in the pitch detector and note extractor steps, 304 and 305 
respectively. 

10 Silence Discriminator 



Figure 4 is a flowchart of the operation of an exemplary silence discriminator 
204 of Figure 2, the silence discriminator isolating the voiced portion in the 
input wavefomi. The first step is to isolate the voiced portions from the 
15 silence portions of digitised hum wavefomi. By preventing the processing of 
silence portions, it Improves perfomiance and reduces computation. A data 
structure is set up using the syntax of the C programming language. 



struct markers{ 
20 int start; 

int interval; 

}; 



where markers is the struct that marks the start and the interval of the voiced 
25 portion. Thus there is an anray of these markers with segjcount members. 



The necessary parameters are initialised to: seg_count=0, can_start='1 and 
count^O, as shown in 401. The parameter can_start is initialised to '1' to 
signal that a new marker is allowed to be created. This is to prevent creating 
30 markers before an interval of voiced portion is registered. It Is followed by 
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process 402 to compute the short-time energy function of the digitised hum 
waveform. The digitised hum waveform is divided into blocks of equal length. 
The short-time energy, E„ , for each block is computed as: 

^" CAL_LENGTH ? N'^M^ " 

5 

where x(m) is the discrete time audio signal, w(m) is a rectangle window 
function and CAL_LENGTH is the length of window and the width of a block 
of hum waveform. 

10 In order to be adaptive to different recording environments, the threshold, 
thres, is computed as the average of the short-time energy and a count 
number is set, 1=0, as shown in 403. The thres is the average short-time 
energy. This is a reference value used to decide whether the signal at a 
particular time is silence or voiced. With the threshold, the short-time energy 

15 of each block Is tested as shown in 404 and 405. In 404, the current short- 
time energy value, energy(i), is tested to detemnine whether its level is greater 
than or equal to 0.9 times the threshold and, at the same time, the 
can_start=1. If the criteria are met, the process proceeds to block 406, where 
the start of the current block is registered as the start of a voiced portion in 

20 406. The position Is calculated as: 

markers[seg_jcount].stait=i*CAL_LENGTH 
where / is the index of the current short-time energy. 

25 

Furthemnore, the can_set is set to '-1' to indicate that the algorithm is 
expecting a silence portion hence another voiced portion cannot be 
registered. If, in step 404, the criteria are not met, the process goes to step 
405, where the cun-ent short-time energy value, energy(i), is tested to 
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determine whether its level is below 0.5*thrBS and, at the same time, the 
can_start=-1 . This is taken to mean that the beginning of a silence portion 
has been reached and, If these criteria are met, ttiis Is registered as an 
Interval In the voiced portion in step 407. The position is calculated as: 

5 

markers[segjcount]/inten/al-i*CALJL£NGTH-markers[segjsount^^^^ 

Following this, the can_start Is set to 'V again to flag that the registration of 
new marker Is allowed and the seg_count is Incremented as shown in 408. 

10 The outputs of steps 406 and 408, together with the output of step 405 if the 
criteria are not met. rejoin in step 409, which asks if all blocks have been 
tested. If the answer is negative, /, the index of the cunrent short-time energy 
is incremented by 1 in step 410 and the process returns to step 404. The 
processes of steps 404 - 410 are repeated until all the values In the short- 

15 time energy function have been tested. 

Gradient Based Segmentation 

The flowchart of exemplary gradient-based segmentation in this invention is 
20 shown In Figure 5A. The humming signal is similar to an amplitude 
modulated (AM) signal where the volume Is modulated by the pitch frequency. 
The pitch signal is not useful for the segmentation algorithm. Thus, the pitch 
frequency is removed to simplify matters. The envelope detector step 501 
removes the pitch frequency. In this way, only information pertaining to the 
25 variation of volume Is left. The differentiator step 502 processes this variation 
to produce a gradient function and removes small gradient values in the 
gradient function. Finally, the note marker extractor step 503 extracts note 
maricers from the threshold gradient function. A typical humming signal with 
three notes hummed is illustrated in Figure 5B. The outputs of envelope 
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detector, differentiator and note markers extractor are illustrated In Figure 5C, 
5D and 6E respectively. 

Envelope Detector 

5 

Figure 6 shows a flowchart for an exemplary envelope detector that Is utilised 
in the gradient-based segmentation as shown in 501. The envelope detector 
consists of two steps: full wave rectification (processes 601 through 605) and 
a moving average low-pass filter. 

10 

The rectifier is simple. In step 601 a count of points in the signal, /, is set to 
"1=0". Following step 602 determines if the signal level at the cun-ent signal 
point is greater than or equal to zero. If it is not, then, in step 603, the 
envelope level for that point is set to the negative of the cun-ent signal level 

15 and i is incremented by 1 in step 605. If the current signal point Is greater 
than or equal to zero, then. In step 604, the envelope level for that point is set 
to the actual signal level and I is incremented by 1 again in step 605. Step 
605 is followed by step 606, which detemiines if "kLEN", where LEN is a 
sample number, chosen here to be 200. If it is, then the process reverts to 

20 step 602. If it is not, then the process goes on to the filter. 

The low pass filter is Implemented by a simple moving average filter to obtain 
a smooth envelope of the discrete time audio signal. In spite of its simplicity, 
the moving average filter is optimal for common tasks such as reducing 
25 random noise while retaining a sharp step response. This property is ideal for 
this invention, as it is desirable to reduce the random-noise-like roughness 
while retaining the gradient. As the name implies, the moving average filter 
operates by averaging a number of points from the discrete signal to produce 
each point in the optimal signal. Thus It can be written as: 
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1 ENVLEN-l 

Where x(0 is the discrete time audio signal with LEN samples. y(t) Is the 
envelope signal of jc(Oand ENVLEN Is the number of points in the average. 
The ENVLEN is chosen to be 200 in this exemplary embodiment. 

5 

The process 607 initialises the necessary parameters "temp", T and '7 to 
zero to start the filtering proper. Before proceeding to filtering, the process 
608 makes sure that the filter operates within the confine of the discrete time 
audio signal, by checking that the sum 7+/<LEAr. The processes 609 and 

10 610 compute the summation of all data after the cunrent value. In particular, 
step 609 provides an updated temporary summation, with "temp=temp+p+jr- 
The average value of the envelope for all T within the sample Is computed as 
shown in 611. ''eny\\]=temp/ENVLEN". Step 612 tests whether the process of 
steps 608 to 611 has been repeated for all data In. the Input buffer and only 

15 when It has does the envelope process end. The T and "j" are incremented 
as show In 609 and 610 respectively. The «++/' Is a pre-increment which 
means j is incremented between testing the condition. 7++" is a post- 
Increment, which means 7" is Incremented after execution of the equation 
shown in steps 610. 

20 

Differentiator 

The flowchart of an exemplary differentiator is shown In Figure 7. The 
difl'erentlator consists of two steps: gradient computation and moving average 
25 low-pass filter. The differentiator processes the envelope produced by the 
envelope detector to generate a gradient function. The algorithm only 
computes the gradient values within the voiced portions marked by the 
markers produced by the silence discriminator. The gradient function 
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essentially describes the changes of the input signal. This can be computed 
by: 

^(f) „ yjt + GRADLEN) - y(t) 
dt GRADLEN 

5 

where y(t) is the envelope signal and GRADLEN is the deviation of f to the 
next point. The GRADLEN Is chosen to be 20 in this exemplary embodiment. 

The process is initialised in step 701 . The index T keeps track of the segment 
10 that is being processed. The index 7" keeps track of the number of points 
within one segment is processed. A decision 702 prevents the overflow of the 
buffer that contains the envelope. "I+Gradlen" is tested against TEN" to 
prevent overflow of the buffer as shown in 702. The gradient is computed by: 

L 

15 where "L" is the step length, for instance 100. Therefore when there is an 
overflow, in step 703 the x(l+L) is set to zero. When there is no buffer 
overflow, the gradient is computed according to this above equation in step 
704. The computation in process 703 caters to the case when the gradient to 
be computed is near the end of the buffer. The step 705 checks whether all 

20 the gradients within the T voice segment are computed. If it is true, it will 
proceed to step 706, else to decision 702. The step 706 increments the 7" to 
process the next voiced segment. The T is initialised to zero to start from the 
beginning of the segment. The decision 707 will check whether all voiced 
segments have been processed. It will proceed to decision 702 if not ail 

25 voiced segments are processed. 

The process 708 initialises the necessary parameters for the filtering 
operation. The filter smoothens the gradient to reduce roughness. The index 



wo 2004/034375 



PCT/SG2003/000232 



18 

of the buffer is tested as shown in 709 to prevent buffer overflow. The moving 
average filter is chosen to smoothen the gradient function. The filter is only 
applied to the voiced portions to reduce computation. The filter length is 
defined as FLEN and all data after the current value are summed as shown in 
710. The index k is tested if it is greater than FLEN as shown in 711. The 
FLEN is chosen to be 200 in this invention. When the FLEN Is reached, the 
gradient, grad, is updated as shown in 712. The process is repeated for all 
points inside the voiced portions, as shown in 713. The processes 709 
through 714 are repeated until all voiced portions are processed. 

Note Makers Extractor 

Ideally, there is only a pair of positive and negative gradient peaks to mark 
the start and end of a note. However, human humming is not ideal and the 
IS problem is further complicated by the presence of expression that causes the 
amplitude in a particular note to change. Thus the note markers extractor has 
to remove invalid gradient peaks based on predefined criteria. These criteria 
are derived from the assumption that each note must be mari<ed by an attack 
and followed immediately by a decay. Anything in between is considered a 
20 false alarm and has to be removed. Figure 8 shows an example that 
illustrates the idea. Figure 8 illustrates exemplary criteria for selection of 
legitimate attack and decay. The criteria of selection of the legitimate attack 
and decay are based on the idea that there is only one attack and decay for 
each note. The 1306 marker is the legitimate attack as it is the first marker 
25 detected. Since decay mariner is expected, the 1307 marker is a false attack. 
Further down, the 1308 marker is temporary considered a decay marker. It 
will be a legitimate decay mari<er if an attack marker follows it. However, a 
decay marker 1309 follows it. Thus, marker 1308 is discarded and the 
marker 1309 is temporarily considered a decay marker. The detection of the 



5 
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attack marker 1310 means that marker 1309 can be formally registered as a 
legitimate decay marker. 

The flowchart in Figure 9 shows an exemplary implementation of the above- 
5 mentioned technique to remove redundant markers. The note markers 
extractor removes redundant ON/OFF markers and registers a set of 
legitimate note markers. A gradient thrashold module 1001 is first called to 
remove small gradient values generated by the differentiator 212. It produces 
a train of ON/OFF pulses. An edge detector function is called to search for 
10 edges from the ON/OFF pulse starting from location 0 as shown in 1002. 
With the location of the nearest marker, the necessary parameters are 
initialised as shown in 1003. In the process 1003, pos and pg are: 



Parameter 


Definition 


pos 


Location of the legitimate attack and decay in the gradient 
array. 


pg 


The gradient value of the legitimate attack and decay. 



The algorithm enters a loop to search and remove all redundant markers as 
shown In 1004 through 1015. The next edge is detected using the edge 
detector starting from the location of the edge found In the last search as 
shown in 1004. The test 1005 ensures that the edge detector has found an 
edge. The 1007 tests for the case when an attack marker is detected while 
an attack marker is registered in the previous iteration. In this case, the 
attack marker detected is discarded and the index is incremented to the 
location of the attack marker as shown in 1011. The 1008 tests for the case 
when a decay marker is detected and an attack marker is detected In the 
previous iteration. Thus, the decay marker detected is registered as a 
legitimate decay marker as shown In 1012. The 1009 tests for the case when 
a decay marker Is detected but a decay marker Is registered at the previous 



15 



20 



25 



wo 2004/034375 



PCT/SG2003/000232 



20 

iteration. Tlius, the current detected marker replaces the previous one as 
shown in 1013. Finally, the 1010 tests for the case when an attack marker is 
detected and a decay marker is detected in the previous Iteration. Therefore, 
the attack marker Is registered as shown in 1014. At a time when the edge 
5 detector Is unable to find any edge, there is a final registration of markers for 
those still pending, as shown in 1006. Since there are no more edges, the 
process 1006 breaks out of the loop and continues to the process 1016. The 
seg_count is calculated as the half of the total number of markers registered, 
as shown in 1016. The processes 1017 and 1018 update the markers struct 
10 with data from pos. 

Gradient Threshold 

Figure 10 shows a flowchart of a simple method to remove the unwanted 
15 small gradient values. The gradient values are tested as shown in 901 . If the 
absolute value Is less than GRADTHRES, it is set to zero as shown in 904. If 
the value is greater than GRADTHRES and positive, it will be set to a positive 
number. If the value is greater than GRADTHRES and negative. It will be set 
to a negative number. Here +10 and -10 are used respectively as an 
20 example. This process Is shown in 902 through 905. In the end, the gradient 
threshold function will produce positive and negative pulses such as those 
shown in 1301 through 1305. 

Edge Detector 

25 

The On/OFF pulses as shown in Figure 8 symbolise the location of high 
gradients. The positive going edges of the pulses as shown by 1301 and 
1302 are the location where gradient values transit from low to high. On the 
other hand, the negative going edges of the pulse as shown by 1301 and 
30 1302 are the locations where the gradient transit from high to low. Thus the 
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negative going edge of the ON pulse is the turning point of the increasing 
envelope to a level value. The negative going edge of the ON pulse is 
detected using the edge detector to obtain the ON markers such as those 
shown in 1306 and 1307. Similarly, the positive going edge of the OFF pulse 
5 is detected using the edge detector to obtain the OFF marker such as those 
shown in 1308 and 1309. 

Figure 1 1 is a flowchart of an exemplary pulse edge detector. The pulse 
edge detector detects the next positive or negative edge starting from the 

10 location specified by start. The process 801 initialises the search index, /, to 
the desired start location. The ps is set to -1 to signal that no previous 
transition is detected. A non-zero gradient and ps=-1 means that this is the 
first time an edge is found as tested in 802. Therefore, ps is set to 1 to signal 
that the first edge is detected as shown in 804. When the gradient value is 

15 zero and ps=1, the second edge is detected as tested In 803. This Is a 
negative going edge for ON pulses and positive edge for OFF pulses. Having 
detected this edge, the current search index will be return as the edge 
detected as shown In 808. The processes from 802 through 805 will repeat 
until all data are exhausted. If the all data are exhausted as tested in 806 and 

20 no edge Is detected, a -1 will be returned as in 807. 

Pitch Detector 

The pitch detector 202 detects the pitch of all note registered in the markers 
25 data structure. Every note interval is divided into blocks that consist of PLEN 
samples. The PLEN is chosen to be 100 in this invention. Thus the pitch 
detection range for an 8KHz sampled audio signal is between 80 to 8KHz. 
The signal in each block is assumed to be stationary and the pitch 
(frequency) is detected by autocorrelation as shown below: 
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1 PLBN-n-l 

where xik) is the discrete time audio signal. 

With this equation, a collection of pitch values that belong to the same note 
5 might be found. In an ideal case, these values are identical. However, the 
autocorrelation pitch detector is sensitive to harmonics that cause en-ors. 
Furthemiore, the hummer might fail to maintain the pitch within the duration of 
a particular note. 

10 Figure 12 shows the flowchart of an exemplary pitch detector. The process 

1101 computes the square of the Input data. The pitch detector is an 
autocorrelation-based pitch detected with some modification. The processes 

1102 through 1114 compute the normalised autocorrelation function and find 
the pitch values of each block in a note. 

15 

A data structure is set up as described below using the syntax of the C 
programming language. 

struct hum_des{ 
20 int pitch; 

int start; 
int interval; 

}; 

25 where markers is the stmct that marks the start and the interval of the voiced 
portion. Thus there is an anray of these marl<ers with notejcount members. 
The position and interval of a note are registered as: 
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hum_des[i].start=marker[i].start 
hum_cleslj].inteval-markeiti].interval 

where j is the index and 0^j< total number of markers. 

5 

The pitch values detected may vary due to the failure of a user to maintain the 
pitch within a single note. The FIndDom function as shown in 1116 finds the 
dominant pitch value, in this invention, the detected pitch values are 
corrected to the nearest MIDI number in 1118. The IVIIDi number is computed 
10 as: 



floor 

hum _ des[j\.pitch = 49 + 



440 ^ 



log2 



The floor(x) function returns a floating-point value representing the largest 
integer that is less than or equal to x. The process Is repeated until all notes 
15 In the input data have their pitch detected as shown in 11 1 9. 



Dominant Pitch Detector 



The function of a dominant pitch detector is to collect statistics from the 
20 collection of pitch values to find tiie prominent pitch values. In this invention, 

the k-mean clustering method is selected to find the prominent pitch values. 

The k-mean clustering method does not require any prior knowledge or 

assumption about the data except for the number of clusters required. 

Determining the number of clustering is problematic in most applications. In 
25 the current invention, the clustering algorithm only needs to cluster the pitch 

values into two groups: the prominent cluster and the outiier cluster. 
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Figure 13: is a flowchart of an exemplary dominant pitch detector (step 1117 
of Figure 12), which uses a l<-mean clustering algorithm that classifies the 
pitches into these two groups. The l«-mean clustering is an iterative algorithm 
for clustering data to reveal the underiying characteristic. The number of 
5 pitches is tested to check if it is greater than 3, as shown In decision 1202. 
The lower and upper 20% of the data are discarded to avoid portions of the 
note that are unstable as shown in 1204. All the pitches will be used for the 
computation if the number of pitches is less than 3. This is attained by setting 
"lower^O" and "upper" to the number of pitches as shown in 1203. The 
10 centres of the two clusters are initialised to the maximum and minimum 
values of the data set as show in 1201 through 1210. The index "j" »s set to 
the lower, as shown in 1205. The process 1211 initialises the necessary 
parameters and saves the cunrent centres for comparison at a later stage. 

15 The pitch values of the note under test are contained in the array pitch. The 
process 1212 compares the absolute distence of the pitch value from the two 
centres. The pitch value is added to the accumulators called, tempi or temp2 
depending on the result of the comparison as shown In 1213 and 1214. This 
process repeats until ail the pitch values in the note are tested as shown in 

20 1215. The new centres are computed and the member counts are 
incremented as shown In 1218 and 1219. They are the average of the 
member pitch values. The processes 1220 and 1221 test if the two centres 
change. If the two centres do not change, the iteration stops immediately. If 
there are changes in any of the centres, the iteration of the processes from 

25 1211 through 1221 repeat unfit the maximum number of loops (MAXLOOP) 
has been reached. The maximum number of loops is 10 in this exemplary 
embodiment. 

If the numbers of members of the two centres is close, as tested in 1223, the 
30 average of the two centres is returned as the dominant pitch. If they are not 
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close enough, the centre with the larger number of members Is retumed as 
the dominant pitch as shown in 1225 through 1227. In this way, the cluster 
with the highest number of members is classified as the prominent cluster 
while the other cluster is classified as the outlier cluster. The pitch of the note 
5 is set to the centre of the prominent cluster. 

It is in fact possible for the invention to work without the silence discriminator. 

Note extraction is a simple module to gather infomiation from note marker 
10 generator and pitch detector. It then filled a structure that describe the begin 
time, duration and the pitch value. Feature extraction converts the note 
descriptors to feature that are used by the search engine. The cun-ent feature 
is the melody contour that is specified in the MPEG-7 standard. The 
description generation is an optional module that converts the feature to a 
1 5 format for storage or transmission. 

Effects of Invention 

The invention achieves the conversion of human (or animal - e.g. dolphin et 
20 al) humming, singing, whistling or other musical noises to musical notes. The 
gradient-based segmentation goes beyond the traditional segmentation 
method that relies on silence. The modified autocorrelation-based pitch 
detector can tolerate a user's failure to maintain pitch within a single note. 
This means that the user can hum naturally without consciously trying to 
25 pause between notes, which may not be easy for some users with little 
musical background. 

While exemplary means of achieving the particular component processes 
have been illustrated, other means achieving similar ends can readily be 
30 incorporated. 



